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Making a system state follow a prescribed trajectory despite fluctuations and errors commonly 
consists in monitoring an observable (temperature, blood-glucose level . . . ) and reacting on its con- 
trollers (heater power, insulin amount . . . ). In the quantum domain, there is a change of paradigm in 
feedback since measurements modify the state of the system, most dramatically when the trajectory 
goes through superpositions of measurement eigenstates. Here, we demonstrate the stabilization of 
an arbitrary trajectory of a superconducting qubit by measurement based feedback. The protocol 
benefits from the long coherence time (T2 > 10 fis) of the 3D transmon qubit, the high efficiency 
(82%) of the phase preserving Josephson amplifier, and fast electronics ensuring less than 500 ns 
delay. At discrete time intervals, the state of the qubit is measured and corrected in case an error is 
detected. For Rabi oscillations, where the discrete measurements occur when the qubit is supposed 
to be in the measurement pointer states, we demonstrate an average fidelity of 85% to the targeted 
trajectory. For Ramsey oscillations, which does not go through pointer states, the average fidelity 
reaches 75%. Incidentally, we demonstrate a fast reset protocol allowing to cool a 3D transmon 
qubit down to 0.6% in the excited state. 



The coupling of a quantum object to an environment 
is essential to enable its observation and manipulation. 
Yet, the mere existence of this coupling induces decoher- 
ence towards pointer states stable under monitoring of 
the environment pQ. There is thus a limiting timescale 
for the faithful preparation of a qubit in an arbitrary 
state, or its control along a given trajectory in Hilbert 
space. As a part of the environment, an observer ex- 
tracts information on the object and contributes to this 
timescale. However, if the observer acquires information 
faster than the uncontrolled part of the environment, it 
is possible to use it through a feedback process and sta- 
bilize permanently a given trajectory or state PH5]. Su- 
perconducting qubits in cavities offer a test bed for these 
concepts, as well as good candidates for practical appli- 
cations j5] |7] . Recently, persistent Rabi oscillations have 
been demonstrated via analog measurement based feed- 
back using continuous weak measurement of a qubit [5], 
and qubit reset via digital measurement based feedback 
using projective measurements has been performed 
In this article, we demonstrate a simple protocol to sta- 
bilize any trajectory of a single qubit using a stroboscopic 
digital feedback based on strong measurement |TD] . Dur- 
ing the manipulation of the qubit, its state is measured in 
a nearly projective manner at specific time intervals and 
a correcting control sequence is triggered conditionally on 
the outcome so as to correct its trajectory from the errors 
due to decoherence and relaxation. The efficiency of the 
trajectory stabilization relies on the rapidity to measure 
and react compared to decoherence. In order to mini- 
mize these timescales, we used a phase preserving quan- 



tum limited amplifier |11H13| and a Field Programmable 
Gate Array (FPGA) able to output a conditioned drive 
pulse less than 500 ns after readout. 

The superconducting qubit follows the design of the 
"3D transmon" developed in Ref. [14J. A single aluminum 
Josephson junction, connected to two antennas of 0.4 mm 
by 1 mm each, on a sapphire substrate, is embedded in an 
empty bulk aluminum cavity whose first coupled modes 
are at lu c /2tt = 7.748 GHz and 13 GHz. External cou- 
pling rates to the first mode Ki n /2ir = 0.34 MHz and 
K out /2ir = 1.49 MHz are chosen of the same order of 
magnitude as the inverse feedback delay (500 ns) and in- 
ternal losses are negligible on these scales. The cavity is 
anchored to a dilution fridge below 30 mK [15]. Spec- 
troscopic measurements give a qubit frequency uj eg — 
lj c — A — 2tt x 3.576 GHz differing from the next tran- 
sition by an anharmonicity (ui eg — uif e )/2ir — 198 MHz. 
The relaxation time T-y = 28 fis corresponds to the Pur- 
cell limit [16 , and pure dephasing time is = 14.5 /Lts 
(Fig#). 



In the dispersive limit u eg — uif e A |17j . the cav- 
ity resonance frequency decreases by 2\ when the qubit 
goes from ground to excited states, and the dispersive 
shift here is xj2-K — 0.78 MHz. The measurement is 
strongest at readout frequency uj r — uj c — x which min- 
imizes the overlap between the two coherent states cor- 
responding to the \g) and \e) qubit states (Fig. [IJ. In 
the experiment, 1.2 [is long square measurement pulses 
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Figure 1: (a) Schematics of the experiment. The state of a 
3D transmon qubit is entangled with the phase of a coherent 
field transmitted through the cavity at frequency uj r = ui c — \- 
It is amplified by a Josephson Parametric Converter (JPC) 
and its complex amplitude a is measured and averaged by 
digital demodulation using an FPGA board (sensing and con- 
trol). The drive at qubit transition frequency cu eg is modu- 
lated by the sum of a predetermined waveform and of a con- 
ditional one generated by the FPGA board (actuation). If 
the measurement points towards state \e), the actuator gen- 
erates a 7T pulse to get the qubit back in \g). (b) Expected 
complex amplitude of the field in the cavity averaged over the 
measurement time T mcas represented as the rod of a lollypop 
in the Fresnel space for both qubit states. The typical devia- 
tion due to vacuum fluctuations of the field in the A r m = 11 
averaged modes is represented by the lollypop radius. The 
limited measurement efficiency (77 = 67%) only slightly in- 
creases the observed deviations by g' 1 ^ 2 — 1 = 22% (purple 
rings). (c,d) Probability density with the JPC OFF (c) and 
ON (d) , extracted from 10 6 measurement outcomes when the 
qubit is prepared in states \g) or |e) with equal probability . 
Each outcome is the complex averaged amplitude of the field 
inside the transmon cavity at uj r . The halved probability den- 
sity corresponding to the preparation of \g) only (resp. |e)) is 
plotted in blue (orange) together with the projections along 
the real and imaginary axes. Turning on the pump of the am- 
plifier (d) results in a great enhancement of the measurement 
fidelity compared to the case without (c). 



are sent through the cavity. The amplitude of the read- 
out field inside the cavity can be calibrated from the 
measurement induced dephasing MS c\ function of read- 
out power leading to 1.4 photons on average [15 . The 
outgoing signal is amplified during these 1.2 /is using a 
Josephson Parametric Converter |llfTT5] with 22 dB of 
gain over 6 MHz (Fig. [I] and (TH]) and following ampli- 
fiers before being down-converted and digitalized using 
the FPGA board input. The board averages numerically 
both quadratures of the signal during the steady part of 
the outgoing pulse only (see Fig. 3a), which corresponds 
to about N m — T mcas (ni n + n out ) = 11 temporal modes 
of 1.4 photons. States \g) and |e) for the qubit lead to 
two almost non overlapping coherent states for the aver- 



age intracavity field \ae~ l9 ) and \ae l9 ) with 9 w 40° as 
expected from tan(0) = 2x/( K in + K out) (Fig. [IJd). With 
an ideal setup measuring both quadratures of the average 
complex field a in the cavity, the variance on a should be 
given by l/y/N m [TS]. In the experiment, the 19% loss 
of signal through the input of the cavity 1 Ktot) and 
the efficiency of the detection setup (82%) degrade the 
signal by only 67% beyond this variance (Fig. [TJ . There- 
fore measuring Im(a) > on the readout field indicates 
a qubit in the excited state |e) with a fidelity beyond 
99.8%, taking aside the expected false counts due to re- 
laxation events during readout. All measurement pulses 
in this article are performed according to this procedure 
and the 0.2% infidelity is neglected throughout. Using 
this setup, it is possible to perform almost projective and 
Quantum Non Demolition (QND) measurements of the 
qubit state much faster than decoherence, a crucial in- 
gredient of measurement based feedback. An illustration 
of the discriminating power of the setup is shown as a 
histogram of measurement outcomes (average complex 
amplitude in the cavity) with the JPC amplifier on or off 
(Fig. Ilp.d) for a qubit starting randomly in state \g) or 
|e). 

Table I: Error in the preparation of \g) using zero, one or two 
resets by feedback when starting in the most entropic state 
or in the thermalized state (effectively at 46 mK). 



reset number 





1 


2 


from(|<,>(fl| + |e>(e|)/2 


50 % 


3.6 % 


1.1 % 


from thermalized state 


2.4 % 


0.7 % 


0.6 % 



As a first demonstration of measurement-based feed- 
back, it is possible to cool down the qubit in its ground 
state, beyond what was done by Riste et al. [3] with a 
phase sensitive amplifier and microprocessor controller. 
Quantum information processing requires such removal 
of entropy during initialization or when correcting for er- 
rors [H], and our method allows to do so without fast 
frequency tuning [20 - 23] , post-selection |241 |2"5] or lim- 
ited coupling rate n < \ [26 . An initial measurement de- 
termines the qubit state. If the outcome points towards 
the excited state (Im(a) > 0), the FPGA controller emits 
a square pulse (Fig. [I]) so as to apply a 7r-pulse around 
Y on the qubit only 500 ns after the first readout pulse 
ends. This reset can be used to increase the repetition 
rate of quantum algorithms or to cool down a qubit. As 
an illustration, the qubit is first prepared in the most 
entropic mixed state p = (\g){g\ + |e)(e|)/2 by either ap- 
plying a 7r-pulse or not (the outcomes are averaged over 
these two possibilities). The probability PL\for the qubit 
to be in state |e) is then measured following zero, one or 
more resets by feedback. 

We found that starting from Pi e \ = 50 %, a single re- 
set brings this level down to Pi e \ = 3.6 % which would 
require to thermalize during 110 (is without feedback. 
Yet events where the qubit relaxes between the middle 
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Figure 2: (a) Evolution of (oz) (dots) and of the coherence | (<jx+i<Jy)\ (circles) when the qubit is prepared in state (\g) + \e))/-v2 
at time 0. The color encodes the time identically in all panels. Lines are exponential fits using coherence time T2 = 11.5 zjs and 
relaxation time T\ — 28 /is. (b) Same evolution represented in the Bloch sphere with a Ramsey frequency ojR y /2n = 100 kHz. 
At each time (color), the outcome of qubit tomography is represented as a dot in the Bloch sphere and in the three orthogonal 
projection planes. The large black circles set the scale of the Bloch sphere extrema. (c) Same evolution as in (a) with 
stroboscopic measurement feedback every 4 /zs. State tomography is only performed outside of the sensing and actuation 
periods. Lines represent results of a simulation without extra fit parameters. When the conditional 7r-pulse occurs (actuation), 
the average purity increases so that the coherences are preserved permanently on average, (d) Same evolution represented in 
the Bloch sphere with a Ramsey frequency LjR y /2n — 10 kHz, instead of 100 kHz for a clearer observation of the trajectory. 
The simulated trajectory is represented as a line only for the time interval chosen in (c) for clarity, (e) Evolution of the 
qubit with the same process as in (c,d) but without actuation. The exponential fit using the same T2 as in (a) indicate an 
average persistent coherence of 18 % without any actuation, (f) In the Bloch sphere, the Ramsey frequency is chosen to be 
uj R y/2-K = 100 kHz. 



of the measurement pulse and feedback pulse limit the 
efficiency of a single reset. Doing a second reset immedi- 
ately after the first brings the qubit much closer to the 
ground state with P^ — 1.1 %. This value does not im- 
prove with additional feedback and is limited mostly by 
the excitation of higher qubit states during the first re- 
set. These higher states are almost empty (0.06 %) when 
starting from a thermalized qubit at P^ = 2.4 % and 
two consecutive resets by feedback cool the qubit further 
down to P| e ) = 0.6 %. These results are summarized in 
Table |T] Note that it is easy to prepare any state with 
similar purity using rotations of the qubit once the qubit 
is in state \g). 

It is also possible to stabilize a state like (\g) + \ e))/y/2 
which is not an eigenstate of the measurement opera- 
tor. First, a 7r/2-pulse is applied to the qubit so as to 
prepare it in (\g) + \e))/y/2 with a drive frequency uj eg . 
At any time t. it is possible to realize the tomography 
of the qubit by averaging the measurement outcomes di- 



rectly ((ctz)), after a fast (64 ns) Rabi 7r/2-pulse around 
Y ((ax)) or after a fast (64 ns) Rabi 7r/2-pulse around X 
((cry)), where (Jx,y,z a re the Pauli matrices. By varying 
the phase of the measurement in time, it is possible to 
represent Ramsey oscillations with any frequency wh„. 
Without measurement based feedback, the Bloch vector 
of the qubit decays exponentially both in Z and in the 
X,Y plane (Fig. [2^i,b). The decay in Z is described by 
timescale T\ = 28 /is while the decay in X, Y is described 
by timescale T2 = 11.5 /is. In order to stabilize persistent 
Ramsey oscillations, a measurement of the qubit is per- 
formed after a 7r/2 rotation every 4 jj,s. The rotation axis 
is chosen so that the measurement outcome should point 
to state \g) in the targeted trajectory and the qubit is 
rotated back to its original state by a —ir/2 pulse. Each 
time the qubit is found to be in the |e) state, the FPGA 
controller performs a fast it pulse (actuation) with a de- 
lay of 500 ns after the measurement ends, which occurs 
after the —ir/2 pulse. Using this stroboscopic measure- 
ment based feedback, Ramsey oscillations are indeed pre- 
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served indefinitely (Fig [2j;,d) . Using optical Bloch equa- 
tions |15j . one can calculate the predicted qubit trajec- 
tory corresponding to this protocol (Fig. |2j;) which is 
consistent with the experiment. Deviations to the exper- 
iment likely originate from the change in measurement 
induced dephasing when the JPC is turned off. The av- 
erage purity Tr(p 2 ) of the density matrix p is calculated 
to be 85% from these simulations, the time averaged 
fidelity F — (V ; targ|p(*)|'^taxg) to the target trajectory 
= (|ff)+e to *»*|e))/\/§is F = 76% and the average 
information quantity 1 — Tr(— plogp) = 0.60 bit. Inter- 
estingly, the sole effect of stroboscopically measuring the 
qubit, without any measurement feedback, induces per- 
sistent Ramsey oscillations, even though with less purity 
(52%), fidelity (56%) and information quantity (0.03 bit) 
(Fig[2j:,f). This is due to the relaxation of the qubit dur- 
ing the measurement period towards state \g) making it 
more probable to reinitiate in state (\g) + \e))/\/2 than 
in state (\g) — \e))/\/2 after the measurement ends. This 
stabilization can be seen as a kind of reservoir engineer- 
ing similar to |27| where the natural qubit decay is used 
as the dissipation source. 

In order to illustrate further the flexibility of strobo- 
scopic projective measurement-based feedback, we have 
also stabilized Rabi oscillations. Although it is possible 
to perform this stabilization using analog feedback on a 
weak, continuous measurement jS], we demonstrate here 
that discrete feedback events are more efficient, mostly 
because it avoids measurement induced dephasing when 
the qubit is far from \g) or |e) [ID]. Without feedback, a 
constant microwave signal at uj eg induces a Rabi oscilla- 
tion of the qubit around cry with decay time Tr = 15.5 ps 
(Fig. [3J3) and frequency set to ujr — 250 kHz. In order 
to make the Rabi oscillations persistent, a measurement 
is performed each time the qubit is supposed to be in 
state \g) (Fig. |3^). The FPGA controller then sends a 
fast correcting 7r-pulse (actuation) at the qubit frequency 
uj eg each time the measurement reveals that the qubit is 
in the excited state. In order to optimize the fidelity of 
the feedback controlled trajectory to the targeted Rabi 
oscillation, the precession angle which is left idle during 
the measurement - Zeno effect freezing the trajectory 
anyway - is briefly accelerated before and after measure- 
ment to compensate exactly for that pause (see Fig. [3^,). 
As can be seen in Fig. [3)3, the Rabi oscillations are indeed 
stabilized permanently with this protocol. Their average 
fidelity to the targeted Rabi oscillation is F — 85%, their 
average purity 80% and their average information quan- 
tity 0.50 bit. The discrete correction events lead to visible 
discontinuities in the trajectories restoring the purity lost 
during the last Rabi period due to decoherence. 

This work illustrates the possibilities offered by 
measurement-based feedback for circuit Quantum Elec- 
trodynamics in the case of a single qubit in cavity. We 
have shown here that fast digital electronics combined 
with efficient detection allow to realize elaborate quan- 
tum control protocols on these systems. Future error 
correction codes will benefit from the malleability of a 
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Figure 3: (a) Pulse sequence for stabilizing Rabi oscillations. 
For a typical period of 4 ps, the lines represents the drive am- 
plitude (green) and expected occupation of the cavity (pur- 
ple). The complex amplitude a of the measurement field is 
recorded only during the steady part of the occupation (red 
area). When Im(a) > 0, a fast ir pulse is applied after a 
total delay of 500 ns (actuation), (b) Black line: decaying 
Rabi oscillation around ov with frequency uir — 250 kHz and 
measured decay time Tr = 15.5 ps. Dots on line: persistent 
Rabi oscillations measured using the pulse sequence described 
in (a), (c) Same measurement as in (b) shown on a smaller 
span for ((Jz) (dots) and (ax) (circles). The targeted Rabi 
trajectory is shown as two lines, dashed for ((Jz) and gray for 



numerical approach where complex filters need to be used 
to protect a quantum algorithm from errors. Any correc- 
tion code needs to evacuate the entropy contained in the 
error, and the reset, presented here, is a valuable tool to 
realize it. Extending these protocols to multi-qubit archi- 
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tectures should enable the preparation and stabilization 
of more complex entangled states and trajectories. 
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